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Abstract 

The missing-mass spectrum obtained in a recoil-free transfer reaction p( Al, 3 He)7r~p'X is 
analyzed. We find that the observed peak structure arises from the coherent contributions from 
two reaction processes in the energy region corresponding to a bound eta (77) meson. In one of 
the processes the intermediate n is captured by the nucleus to form the 77-mesic nucleus 25 Mg i) . 
In the other process, the n does not form 77-nucleus bound state. The interference between these 
two processes has caused the peak of the spectrum to appear at an n binding energy stronger 
than the actual one. Our analysis also indicates that the data are consistent with an attractive 
N*(1535)-nucleus interaction at energies below the 77N threshold. 

Keywords: Nuclear interaction of N*(1535), Eta-mesic nucleus Mg^. 
PACS Numbers: 24.10.Eq, 24.10.-i, 21.10.Dr 



* e-mail address: haider@fordham.edu 
t e-mail address: liu@lanl.gov 



1 



I. INTRODUCTION 



The possibility of an eta (rj) meson being bound in a nucleus by strong-interaction 
force, leading to the formation of an exotic nucleus called eta-mesic nucleus, was first pro- 
posed in 1986 l|. The nuclear binding of an 77 meson is caused by an attractive eta-nucleon 
(7/N) interaction in the threshold region. This attractive 77N interaction is due to the N*(1535) 
resonance being situated not too far above the threshold of the 77 N channel (1488 MeV). It 
was also shown in Ref. [lj] that a bound state is possible if there is a sufficient number of 
nucleons (10 or more) in a nucleus. 

The existence of ?7-mesic nuclei will certainly create new premises for studying 77 
meson and baryon resonances inside nuclei. Extensive interest in exploring this new venue 
is evidenced by the large amount of theoretical work that has been done during the last two 
decades 2|- 16J. Irrespective of the models and formalisms used, from a theoretical point of 
view there is unanimity about the existence of 77-mesic nucleus. 

n 

Recently, the COSY-GEM collaboration [17J searched for 77-mesic nucleus by means 
of a recoil-free transfer reaction p( 27 Al, 3 He)7r~p'X. The kinematics was so chosen that the 
77 produced in the intermediate state is nearly at rest, favoring its capture by the residual 
nucleus 25 Mg. Because of energy conservation, the bound 77 cannot reappear as an observable 
particle in the decay products of the mesic nucleus 25 Mg^. Instead, it interacts, for example, 
with a target neutron resulting in the emission of a nearly "back-to-back" ir~p pair in the 
laboratory. We denote this multi-step reaction as Process M (M for mesic-nucleus) : 

p + 27 A\ -> 77 + 25 Mg + 3 He 

S v ' 

I 

25 Mg„ 

; 

r] + 25 Mg^ (71- +p)+X . 
In order to reduce a large number of background events arising from particles being 
emitted during nuclear cascade process, the COSY-GEM collaboration implemented a triple- 
coincidence measurement among 3 He and the "back-to-back" ir~p pair having the kinetic 
energy spectra of the pion and proton peaked, respectively, at about 100 and 320 MeV fis| . 
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Because of this background reduction, a peak in the missing-mass spectrum in the energy 
region corresponding to bound 77 has been made evident. Upon fitting the spectrum with 
the sum of a background term and a Gaussian (\fb\ 2 + |/ 9 | 2 ), COSY-GEM collaboration 
determined that the peak has its centroid situated at binding energy (—13.13 ± 1.64) MeV 
with a FWHM of (10.22 ± 2.98) MeV [or a half- width T/2 ~ (5.1 ± 1.5) MeV]. 

By performing a bound-state calculation based on scattering length and using on- 
shell kinematics, we find that the above binding energy and half-width correspond to an 
effective s-wave 77N scattering length ao — (0.292 + 0.077i) fm. An exceptional feature of 
this scattering length is its imaginary-to-real part ratio R = Xm(ao) /TZe(a ) = 0.26 only, 
while most of the published theoretical models give scattering lengths (see table |T] of Ref.[19|) 
having R 3> 0.35. In other words, the value of R given by the theories is higher than the 
fitted value by at least 35%. The need to understand the huge difference between theory 
and experiment has motivated the present study. 

More specifically, we will reanalyze the experiment and infer from our analysis the 
nature of the observed peak structure and the qualitative feature of the N*(1535)-nucleus 
interaction. In section II we outline the mesic-nucleus theory to be used in the analysis. 
Detailed analysis is given in section III, and our findings are summarized in section IV. 



II. OUTLINE OF THE MESIC-NUCLEUS THEORY 

The eigenvalue equation of the bound-state of an rj meson in a nucleus is (H +V) \ip) = 
£\ip). The ^-nucleus potential V is complex because the rj — » it channels are open. Hence, 
the eigenenergy £ is also complex and can be written as £ = E bd —iT/2, where E bd (< 0) and 
T are the binding energy and width of the bound state, respectively. The momentum-space 



matrix elements of the leading-order potential are given by 19|]- 20] 



(k>\V\k) = {K'\t vN (W)\K)F(k'-k), (1) 

— * — * 

where k and k' denote the initial and final r\ momenta in the 7/-nucleus cm. frame and 
F(k' — k) is the nuclear form factor. The t V N is the operator for the scattering of rj from 
the nucleon. The variables k, k', and W are, respectively, the initial and final r/N relative 



momenta, and the total energy of the 77N system in its cm. frame. 

Without loss of generality, we use the coupled-channel isobar (CCI) model of Bhalerao 
and Liu jiZlJ to calculate the potential given by Eq.([T]). The reason for this is two-fold. First, 
we have at our disposal the detailed energy dependence of the model which reproduces 
remarkably the observable (7rN Sn phase shifts) in the entire energy region where the nuclear 
binding of an 77 could take place. Second, it was this model that was used to predict the 
existence of 77-nucleus bound states. Furthermore, as has been noted in the previous section, 
the discrepancy between the theoretical and experimental value of R exists for all published 
models, including the CCI model. Hence, we believe the general features of our findings are 
not limited to the model used. 

In the CCI model of Ref . (2J , 



(K>\t vN (W)\K) = Kj2v e (K',A e )A e (W)ve(K,Ae), (2) 

1 

where K is a kinematic factor and I = 0,1,2 are, respectively, the s-, p-, and c?-wave 
77N interactions. The V£ are off-shell form factors of range A^, and Ae are the energy- 
dependent amplitudes. For bound-state problems, p- and d-wave interactions have negligible 
contributions. Consequently, we will only consider the s-wave 77N interaction and omit the 
subscript I. The amplitude A is given by 

Conversely, if A is the amplitude that gives the measured E bd and T/2, then Eq.()3]) can be 
used to derive the energy dependence of the denominator, namely, 

D(W) = (4) 
V ; 2WA 

In the above equations g is the r/NN* coupling constant. Here, N* is the s-wave isobar 
N*(1535) which has a mass between 1525 and 1545 MeV and a Breit-Wigner width from 125 
to 175 MeV [221 ] . For mesic-nucleus calculations, 

W = m v + m N + (B N ), (5) 
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where (B N ) < is the average binding energy of the nucleon. In Eq.(j3]) 

D(W) — W — M N *(W), (6) 

and 

M N .(W) = M° + r(W) + TZe[Tr ed {W)} + i Xm[£ med (H/)], (7) 

with J] med = j]f ree _)_ j] abs . The J]f ree is the N* self-energy arising from its decays to the 77N, 
7rN, and ttttN channels in free space. If , y^ ree (W) denotes its total free-space decay width, 
then \rf ree {W) = -Im[^ ree {W)\. The bare mass M°, coupling constant g, range param- 
eter A needed for the calculation of £^ ree were determined from fitting the experimental 
7rN Sh phase shifts 21[. The Y> abs is the N* self-energy arising from true absorption (or 



annihilation) by nucleons of the pions coming from N* — > 7rN, and tcttN decays. To our 
knowledge, microscopic absorption model that can fit systematically all experimental data is 
still not available. In the literature, Tm\T, abs ) has been estimated in the framework of local- 



density approximation, and 7?.e[E a6s ] is treated as a parameter 23]. It was found that total 
-Im[S a6s ] ~ 35 MeV at W = 1535 MeV and at nuclear density p = p = 0.17 fm" 3 . We 
extend this result to the subthreshold region by using — Zm[E ] = 35(q/y) 3 (p/ p ) 2 MeV. 
Here, y is the 7rN relative momentum at W = 1535 MeV and q the corresponding one at 
W < (m N + m^). The exponent of q is based on the local-density result [23(] on the momen- 
tum dependence of T, abs while the exponent of p is based on the fact that in a nucleus, pion 
absorption involves at least two nucleons. In principle, any subthreshold N*-nucleus inter- 
actions can contribute to the real part of the N* self-energy. We denote these contributions 
collectively as r(W) in Eq.(J7]). 

We emphasize that Xm[E abs ] has the same sign as Im[S^ ree ], as required by the uni- 
tarity of an optical potential. Hence, |Xm[S med ] | > |Xm[S^ ree ] |. The dependence of £^ ree 
and the maximal £ med on W are shown, respectively, as the dash-dotted and solid curves 
in figiU At any given W, a physically meaningful ^4(iy) must always yield an Xm[S med ] 
situated in the "physical zone" bordered by these two curves, which we term the unitarity 
requirement. In fig{T] the left and right vertical lines show the positions of W correspond- 
ing, respectively, to nucleon binding energies (B N ) = —100 and —30 MeV. As one can see, 
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FIG. 1: The energy dependence of — Im[S^ ree ] (dash-dotted curve) and 
+ Im[T l abs ]p = p o ^ (solid curve). The left and right vertical dotted lines indicate, respectively, the 
W corresponding to {B N ) = —100 and —30 MeV. The dashed curve gives the energy dependence 
of —Im[T> med ] in case the experimental spectrum [17] is due solely to mesic-nucleus formation (see 
the text in Section IIIip . 

between these two (B N )'s the lower boundary of Xm[E med ] is nearly constant. 

The quantity r + 7?.e[S a6s ] in Eq.(|7]) represents the real part of the N*-nucleus inter- 
action which we denote as Vn* . Equation (J7|) can then be written as 

M N *(W) = M° + V N *(W) + TZe[J: free {W)] + i lm[E med {W)] . (8) 

Using Eq.® in Eq.(|6]), we obtain 

D{W) = W - ( M° + V N *(W) + Ke[E free (W)] + i lm[^ med {W)} ) . (9) 
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Upon equating the real and imaginary parts of Eq.(j9]), one obtains 

V N *(W) = W-M(W)-Ke[D(W)], (10) 

and 

lm[Y, med {W)} = -Xm[D{W)\. (11) 

In Eq.fllDjl, M(W) = M° + TZe[J: free (W)} and W is given by Eq.([5J. We recall that (B N ) 
is the average nucleon binding energy. Because the nucleons are bound, (B N ) < 0. On the 
other hand, Vjy« can be either negative or positive, depending on whether the interaction is 
attractive or repulsive. 

III. ANALYSIS AND DISCUSSION 

If the centroid of the experimental peak at —13 MeV and half-width of 5 MeV are 
due solely to the formation of the //-mesic nucleus 25 Mg^ via Process M, then an amplitude 
A = —(0.0521 + 0.0099i) fm 2 is required to reproduce the above data. Upon using Eqs.(j3J 
and (TlTj) to solve for Tm[Y, med (W)}, we obtain the dashed curve in figiU As one can see, 
this curve intersects the physical zone at W c± 1125 MeV, which, by Eq.((5]), corresponds 
to a (B N ) ~ —360 MeV. This is clearly an unrealistic value which we regard as a strong 
indication that Process M alone is insufficient in describing the observed spectrum. 

Indeed, the rj produced in the intermediate state can also be scattered by the residual 
nucleus and emerge as a pion, without being first captured by the nucleus. We denote this 
multi-step reaction as Process S (S for scattering): 

p + 27 Al ->• r) + 25 Mg + 3 He 

, " > 

r/ + 25 Mg-+ (tt- +p)+X . 
The essential portion of the reaction dynamics that differentiates the S and M pro- 
cesses, as indicated by upper and lower braces in the corresponding reaction equations, are 
illustrated in figfJl We emphasize that because these two reaction paths lead to the same 
measured final state, they cannot be distinguished by the experiment. Consequently, in the- 
oretical analysis one must take coherent summation of the two amplitudes to account for the 
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FIG. 2: (a) Reaction diagram of f s . (b) Reaction diagram of f M . The wavy and multiple lines 
represent, respectively, the r\ and 25 Mg. The open oval denotes the ^-nucleus interaction V. The 
filled line in (b) denotes the mesic nucleus. 

quantum interference between them. We, therefore, fit the experimental spectrum by using 
the sum of two amplitudes: 

2 



a I fa + A 



a 



< k'\V(E)\k>^ 



< k'\V(E)\ip >< V\V(E)\k > 



'12) 



E - (E bd - iT/2) 

where V is given by Eq.([I]), E = W—m v —m N — (B N ), ijj is the wave function of bound rj, and 
is its adjoint 24j. We have noted that in the threshold and subthreshold regions, rj- nucleus 
interaction is isotropic and that the matrix elements < k'\V(E)\k > are nearly constant for 
k and k' between and 100 MeV/c. Because of these aspects of the 7/-nucleus interaction 
and the experimental selection of events corresponding to rj being produced nearly at rest, 
Eq.( !T2|) can be evaluated at \k\ = \k'\ ~ 0. 

One should note that in Eq. (II 2j) . there is only one parameter a and its role is to just 
adjust the overall magnitude. We emphasize that we used the same V in calculating f s and 
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f M , and that the values of Em and T/2 obtained with this V were kept fixed during the fit. 
Furthermore, we square the sum of the amplitudes, in marked contrast to using the sum of 
the squared individual amplitudes. Hence, interference effects between the amplitudes are 
present in our analysis while they are absent in the COSY-GEM fit. 

Our main goal is to find answer to the following questions: (l)Would it be possible 
to explain the spectrum by having \Eu\ <^Cl3 MeV in f M l (2) Is Vn* attractive or repulsive? 

In our previous work [if]], we predicted the existence of 77-mesic nuclei with < B N >= 
—30 MeV, but without the inclusion of N*-nucleus interaction. In other words, we have set 
both V^* and Tm\Yf lhs \ equal to zero. The use of average nucleon binding energy < >= 



—30 MeV is based on the findings from various studies of meson-nucleus scattering 25|, |26j. 
Applying the same approach to 25 Mg, we obtained £ = E bd - iT/2 = -(6.5 + Hi) MeV. 
Upon using this £ and the corresponding V to calculate f M and f s , we obtained the result 
given in row (a) of table fl] and the spectrum is displayed as the dotted curve in figj3j As one 
can see, the peak position of this curve appears at —10.5 MeV. This 4.0-MeV downward shift 
from —6.5 MeV indicates clearly the importance of interference effect. With the use of V 
obtained with V N * = -24 MeV and -Xm\L abs \ = 0.65 MeV, we obtained E bd = -8.0 MeV 
and T/2 = 9.8 MeV. The corresponding spectrum given by Eq. ffl2l is shown as solid curve in 
figj3j The peak position of this curve is at —12.5 MeV, which agrees with the COSY-GEM 
result of (—13.13 ± 1.64) MeV. The quantitative aspect of this latter calculation is given in 
row (b) of table fl] 

It is interesting to note from table Hthat fits (a) and (b) give same value to the overall 
conversion parameter a. As indicated in section II, the //-nucleus interaction strength is given 
by D(W) and its values are shown in the 6th and 7th columns of the table. The comparison 
of the 7th and 8th columns of the table shows that Im[D](= -Im[S med ]) > — Xm[S^ ree ]. 
Hence, they satisfy the unitarity requirement. We also did calculations using larger Xm[D]. 
However, they led to larger calculated half-widths which worsened the fit. Consequently, 
we conclude that the present data require Xm[S med ] being close to the lower boundary of 
the physical zone of the self-energy (the dashed-dotted curve in fig HI). This in turn implies 



that in the recoilless transfer reaction of Ref. 



17|, Xm[S a6s ] (=lm[E med } -Zm[£ /ree ]) is 
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FIG. 3: Spectra obtained with (a) potential V giving Em — iT/2 = —(6.5 + 7.U) MeV (dashed 
curve), and (b) potential V giving Ebd — iT/2 = —(8.0 + 9.8i) MeV (solid curve). The data are 
from Ref.Q. 

TABLE I: Quantitative details (a is in counts/fm 2 ; all other quantities are in MeV). 



Fit 


a 


Ebd 


T/2 


<B N > 


Ke[D] 


Xm[D] 


-Im[E^ ee ] 


M 


Vn* 


(a) 


4.2 


-6.5 


7.1 


-30 


-164 


68.3 


68.3 


1622 





(b) 


4.2 


-8.0 


9.6 


-30 


-140 


69.0 


68.3 


1622 


-24 



very small, and pion absorbtion takes place mainly in the nuclear surface region where the 
condition p/po *C 1 leads to a very small Im | E a&s | . 
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IV. CONCLUSIONS 



Our analysis shows that two reaction processes are contributing to the observed spec- 
trum of the bound 77 in 25 Mg. The interference between these processes causes the centroid 
of the observed spectrum to appear at an energy stronger than the actual binding energy of 
the rj meson. The effective r]N scattering lengths that reproduced Em and Y/2 used in f M 
associated with the dashed and solid curves of fig]3]are, respectively, (0.226 + 0.094z) fm and 
(0.250 + 0.123?) fm. The corresponding imaginary-to-real part ratios are R=0A2 and 0.49, 
consistent with theories (see Section [I]). We, therefore, explained the apparent discrepancy 
between theory and experiment. 

The present analysis also indicates that the real part of the interaction between N* and 
25 Mg is attractive at energies below the rjN threshold. This latter new nuclear information 
should be of value to nuclear physics studies involving the baryon resonance N*(1535) in 
medium-mass nucleus, such as 25 Mg. 

We emphasize that the existence of S and M processes and the interference between 
them are of a general nature. Consequently, our finding on the effects arising from this 
aspect of reaction dynamics is model- independent. On the other hand, while the specific 
value of Vat. may be mo del- dependent, its sign (or the attractive nature of the interaction) 
is model-independent. This is because the negative sign is required to provide more binding 
which, when combined with interference effects, can lead to sufficient downward shift of the 
peak in the binding-energy spectrum. We invite other researchers having at their disposal 
the detailed off-shell properties of their models to analyze the COSY-GEM data to further 
pin down the value of V^*. 
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